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ABSTRACT 

PSR J0751+1807 is a millisecond pulsar in a circular 6hr binary system with a helium white dwarf 
secondary. Through high precision pulse timing measurements with the Arecibo and Effelsberg radio 
telescopes, we have detected the decay of its orbit due to emission of gravitational radiation. This is 
the first detection of the relativistic orbital decay of a low-mass, circular binary pulsar system. The 
measured rate of change in orbital period, corrected for acceleration biases, is Pb — (— 6.4±0.9) x 10"^"*. 
Interpreted in the context of general relativity, and combined with measurement of Shapiro delay, it 
implies a pulsar mass of 2.1 ± 0.2 Mq, the most massive pulsar measured. This adds to the emerging 
trend toward relatively high neutron star masses in neutron star-white dwarf binaries. Additionally, 
there is some evidence for an inverse correlation between pulsar mass and orbital period in these 
systems. We consider alternatives to the general relativistic analysis of the data, and we use the 
pulsar timing data to place limits on violations of the strong equivalence principle. 
Subject headings: gravitation — binaries: general — pulsars: individual (PSR 0751+1807) 



1. INTRODUCTION 

A fundamental prediction of general relativity is the 
emission of gravitational radiation by binary systems. 
The consequent loss of energy and angular momentum 
from the se systems causes their orbital periods to 
decrease ()Petersiri964j) . Previous to the present work, 
relativistic orbit decay had been detected in five binary 
pulsar systems. Four of these systems consist of a 
mildly recycled pu lsar bound to a second neutron star in 
an eccentric orb it (Tavlor & Wcisbcre 1989: Stair s et all 
l2002HDeich~Kulkarni , 1996: .Kramer et al. 200a . The 
fifth system, that of PSR J1141— 6545, contains a young 
pulsar and a ^IM ^:^ white dwarf in an eccentric orbit 
ilBailes et a,] 1l200l. 

In this paper, we report the measurement of the 
relativistic decay of P SR J0751+1807, a mi llisecond 
pulsar in a 6hr orbit ijLundgren et alJ [l995t) . This 
system differs in several ways from the other binaries 
in which relativistic decay has been detected. First, 
its rotation period of 3.4 ms is an order of magnitude 
smaller than the rotation periods of the pulsars in the 
other systems, implying greater accumulation of angular 



momentum during its binary accretion phase. Second, 
its magnetic field is substantially lower than those of 
the other systems, presumably b ecause of diffe rences 
in orbital evolution and accretion l)Bhattacharval 12002') . 
Third, its orbit is extremely circular, with eccentricity 
under 2 x 10^^, a consequence of tida l circularizatio n 
during the late stages of the secondary (Phin nevlll992(l . 
Fourth, the secondary to PSR J0751+1807 is a low 
mass helium white dwarf, several times lighter than the 
secondary stars in the other systems. 

Measurements of relativistic orbital phenomena in 
binary systems yield constraints on the masses of 
the stellar components of the systems. Neutron star 
mass measurements, in turn, serve as probes of the 
properties of nucle ar ma tter at high densities (e.g., 
Lattimcr fc PrakashI l20?)l . The maximum achievable 
neutron star mass depends on the equation of state of 
nuclear matter. Soft equations of state, expected if the 
core of the neutron star is composed of non-nucleonic 
matter, predict the maximum neutron star mass <2M0, 
while stiffer equations of state allow higher values. 

Reviewing the field several years ago, 
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TABLE 1 
Summary of Observations 



Observatory System Dates Frequency Bandwidth Number Typical RMS 

(MHz) (MHz) of Integration Residual 

TOAs (min) {^ls) 



Arecibo Mark HI 1993.8-1994.4 430 8 1190 0.5 25 

Mark IV 1998.9-2004.1 430 10 1007 3.2 5 

Mark IV 1997.9-2001.5 1410 10 109 3.2 7 

WAPP 2004.1-2004.1 1145 50 562 0.5 10 

WAPP 2004.1-2004.1 1195 50 559 0.5 8 

WAPP 2004.1-2004.1 1395 50 559 0.5 8 

WAPP 2004.1-2004.1 1445 50 559 0.5 8 

Effelsberg EBPP 1997.0-2004.6 1409 56 490 7.0 7 



iThorsett fc Chakrabartvl l)1999(l found all radio pulsar 
mass measurements to be consistent with a narrow 
underlying Gaussian mass distribution, 1.35 ± 0.04 M©. 
More recent observations have found (1) neutron star 
masses in neutron star-neutron star binary systems lie 
in the range 1.18 to 1.44M0; (2) PSR J1141-6454, a 
young pulsar in an eccentric orbit with a white dwarf 
secondary, has mass I.SOMq; and (3) neutron stars 
in circular neutron star-wide dwarf systems may have 
some what larger mass es, but the uncertainties are large 
fe.g.. lNice eTal]l2005l) . 

Neutron star masses can also be measured in X-ray 
binary systems. Spectroscopic and eclipse observations 
typically yield neutron star mass measurements in the 
range 1.1 to 1.5 M0, but there is evidence for a few more 
massive neutron star s: 1.78 ± 0.2 for the neutron 
star in Cygnus X-2 l|Orosz k Kuulk crs 199 3), 2.44 ib 
0.27 M for the neutron star in 4U1700-37 ( Clark et all 
I2002|) a nd 2.27 ± 0.17 for the neutron star in Vela 
X-1 dQ uaintrcU ct al. 2003). Se e the review papers by 
ILattimer fc Prakash (,2004) and ICharles fc Coi l|2005fl 
for further references. 

We undertook high precision pulse timing observations 
of PSR J0751+1807 with the goal of measuring its orbital 
decay, both to test the expected relativistic behavior of 
its binary system and to deduce the pulsar and secondary 
star masses. In f]21we describe the observations in detail. 
In fJ31 we present results of fitting pulse timing models 
to the data. In 21 use our measurements to place 
constraints on theories of relativistic gravity. In ^ we 
discuss the pulsar mass and the evolution of the binary 
system. In ^ we summarize the major points of the 
paper. 

2. OBSERVATIONS 

We observed PSR J0751-I-1807 for more than ten years 
using the 305 m telescope at Arecibo and the 100 m 
telescope at Effelsberg. The bulk of the Arecibo data 
were collected in campaigns of duration ^^1 week over 
which all orbital phases were observed. Such campaigns 
were undertaken in November 1993, May 1999, May 
2000, June/July 2001, May 2003, and January 2004. 
Additional data were collected at some other epochs 
in 1993-1994 and 1999-2000 (Figure [TJ. Effelsberg 
data were less concentrated, with observations spread 
evenly over several years and a single major campaign 
in February 2000. The two types of data complement 



each other, with campaign data especially useful for 
measuring orbital elements, and with less concentrated 
observations important for determining pulse spin-down 
and astrometric parameters. Ultimately, all 5035 pulse 
arrival times from both observatories were combined and 
simultaneously fit to pulsar timing models O- 

2.1. Arecibo Observations 

Observations at A r ecibo in 1993 and 1994 are described 
in ILundgren et alJ ljl995() . The Princeton Mark III 
data acquisition system collected signals across a 8 MHz 
passband at 430 MHz using a 32-channel filter bank 
spectrometer with 100 fi time constants. Observations 
from 1997 to 2004 employed the Princeton M ark IV 
coherent dedispersion system llStairs et al.ll2000() . Most 
observations were made at 430 MHz with a 5 MHz 
passband, but some were made at 1410 MHz using a 
10 MHz passband. Observations typically lasted 29 
minutes and were analyzed in blocks of 190 s. 

The January 2004 Arecibo campaign included 
extensive observations at radio frequencies 1120 to 
1470 MHz. These data were collected with four 
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Fig. 1. 



Residual pulse arrival times after removing the best 
fitting pulse timing model. Points are segregated by observatory, 
frequency, and data acquisition system. Note the difference in 
vertical scales. 
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TABLE 2 
Timing Model Parameters^ 



Basic Timing Model (Three Post-Keplerian Parameters) 



Ecliptic longitude, A 116?33362028{2) 

Ecliptic latitude, /3 -2?807548{2) 

Proper motion in A, fi\ = cos l3{d\/dt) (mas/yr) —0.35(3) 

Proper motion in /3, /ip (mas/yr) ~6(2) 

Parallax (mas), n 1.6(8) 

Rotation frequency, uq (s'^) 287.457858630106(2) 

Rotation frequency derivative, ui (s'^)*^ -6.4337(4) X 10"^'' 

Epoch, to (MJD [TDB]) 51800.0 

Dispersion measure, DMo (pccm-3)': 30.2489(3) 

Dispersion measure derivative, DMi (pccm~^yr^^) ... —0.00017(1) 

Orbital period, Pt (days)'' 0.263144266723(5) 

Projected semi-major axis, x (It-s) 0.3966127(6) 

Eccentricity, e 5(11) X 10~^ 

Time of ascending node, tasc (MJD [TDB]) 51800.21573411(2) 

Orbital period derivative, P5 (unitless) -6.2(8) X lO"'* 

Shapiro parameters, r and s see Figure 01 



General Relativistic Timing Model (Two Post-Keplerian Parameters) 



Cosine of inclination angle, cosi OAl_Qg^ 

Pulsar mass, mi (Mq) 2.1(2) 

Secondary mass, m2 (M0) 0.191(15) 



''Figures in parentheses are 68% confidence uncertainties in the last digit quoted. 

Observed value, not corrected for acceleration biases; see Table 13 
'^Formal uncertainty in the timing fit. No attempts v^ere made to correct for pulse shape evolution over 
frequency. 



Wideband Arecibo Pulsar Processors (WAPPs). Each 
WAPP calculated autocorrelations of 192 lags across 
a 50 MHz passband in each sense of polarization. 
Autocorrelations were accumulated for 32 ^s and then 
summed into folded pulse profiles with 256 bins across 
the pulse period. The four WAPP passbands were 
centered at 1145, 1195, 1395, and 1445 MHz. Data from 
each of the four passbands were reduced independently. 

Each data acquisition system folded the signals modulo 
the pulse period using a precomputed ephemeris. A pulse 
time of arrival (TOA) was calculated from each data 
block by fitting a high quality template profile to the 
data profile and adding the phase offsets thus measured 
to the observation start times. A suitable number of 
pulse periods was added to the TOA to make it fall near 
the center of the observation. Start times were referenced 
to the observatory time standard and later corrected to 
UTC and, ultimately, TT(BIPM03). 

2.2. Effelsberg Observations 

At Effelsberg, timing data of PSR J0751-M807 have 
been collected since 1997, with observations made 
approximately once a month using a 1300—1700 MHz 
tunable HEMT receiver at a center frequency 
of 1409 MHz. All observations employed the 
Effelsberg-Berkeley Pulsar Processor (EBPP) as 
the data acquisition system, correcting the dispersion 
smearing of the signal using coherent dedispersion 
l|Ra,nker et a.1 .111991 . 

In total power mode, the EBPP provided 32 channels 
across a total of 56 MHz for each sense of circular 
polarization. The output signals of each channel were 
fed into digital dedisperser boards for coherent on-line 
dedispersion and were synchronously folded at the pulse 



period over 7 niin integration times. The signal-to-noise 
ratio varied between 3 and 10 for the individual 
integrations, depending on interstellar scintillation. As 
with the Arecibo data, TOAs were calculated by fitting 
the data to high quality template profiles and calculating 
the phase offset relative to the observation start times. 
Start times were referenced to the observatory hydrogen 
maser clock, corrected to UTC(NIST) using the signals 
from the Global Positioning System (GPS), and finally 
corrected to TT(BIPM03). 

3. TIMING ANALYSIS 

3.1. Basic Timing Model (Three Post-Keplerian 
Parameters ) 

The TOA measurements were fit to pulse timing 
models using the tempo software package^. The basic 
timing model had 23 parameters and 5012 degrees of 
freedom. The results of the fit is summarized in Table|21 
and residual arrival times after removing the timing 
model are plotted in Figure ^ Pulsar rotation was 
parameterized by spin frequency, uq, and its derivative, 
vi. Earth motion was m odeled by the JPL DE405 
Ephemeris l|Standishl IT991 . To minimize covariances 
between components of position and proper motion, 
astrometric results are presented in ecliptic coordinates, 
calculated by rotating the coordinates of the solar system 
ephemeris about the x axis by the obliquity of the Earth, 
eo = 84381'.'412. Arbitrary time offsets were fit between 
data sets collected at different observatories, at different 
frequencies, or with different equipment. Interstellar 
dispersion was modeled with a dispersion measure, DM, 
changing linearly with time, DM(t) = DMq + DMit, 

^ http:/ /pulsar. princcton.cdu/tempo 



4 



while dispersion within the solar system was calculated 
using an analytical model of the solar wind f >l3.5|l . 

Orbital motion was calculated using the "ellI" 
relati vistic orbital model (Appendix A of iLange et al.l 
1200 A combination of five Keplerian and three 
post-Keplerian orbital elements were fit to the data. 
The Keplerian orbital elements were the orbital period, 
Pt,; the semi-major axis projected into the line of sight, 
X — (aisini)/c, where i is the inclination and c is the 
speed of light; the time of passage through the ascending 
node, tasc! and two Laplace-Lagrange parameters, 77 and 
K, which parameterize the eccentricity. Since neither 
of these were statistically significant, we report total 
eccentricity, e = (ry^ + k^)^/^; the value in Table |21 should 
be taken as an upper limit. 

The three post-Keplerian parameters were the time 
rate of change of orbital period, Pb, and the shape and 
range parameters of Shapiro delay, s and r. The latter 
are defined by the perturbation in arrival time, Ai, 
imposed by Shapiro delay; for a nearly circular orbit, 
this is 



-2rln{l - ssin[(27r/Pb)(i - tasc)]}, (1) 



where t is the pulse arrival time. 

The best fit parameters are given in Table El under the 
heading "Basic Timing Model." The most important 
new measurement is the relativistic period decay rate. 



Ph 



^(6.2 ±0.8) X 10 



-14 



(2) 



discussed further below. The two Shapiro parameters, r 
and s, are highly covariant, and the confidence regions 
are non-elliptical. The values allowed by the timing data 
are shown in Figure El 

3.2. The Basic Timing Model under General Relativity 

The basic timing model provides an excellent fit 
to the data without reference to a specific theory of 
gravitation. The implications of this for gravitation 
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Fig. 2. — Constraints on Shapiro delay parameters r and s from 
the basic timing analysis. Inner and outer contours delimit 68% 
and 95% confidence regions. 



theory are discussed below ([2J- Here, these results are 
interpreted in terms of general relativity. 

Under general relativity, the post-Keplerian 
parameters are related to the masses of the pulsar, 
mi, and the secondary star, 7712, and the inclination of 
the orbit: 



1- 



73 



37 



(1- 



and 



\7/2 ^0 



T0TO2, 



s = smi, 



24 96 
mi m2 

(mi +m2)-'-/'^ 



(3) 



(4) 



(5) 

,3 _ 



where masses are in solar units and Tq — GMq/c 
4.925490947 x 10^^ s. The masses and inchnation are 
related by the Keplerian mass function. 



_ (m2sini)^ 
Ji = 



(mi -I- m2)^ 



X 



(6) 



Under general relativity, the measured values of orbital 
decay and Shapiro delay yield the constraints on mi, 
m2 , and i shown in Figure Evidently, meaningful 
constraints on the individual masses can be attained 
only by combining the Shapiro delay and orbit decay 
measurements. Given this circumstance, it is useful to 
directly incorporate general relativity into the timing 
analysis and to explicitly consider mi and m2 as 
independent variables. 

3.3. General Relativistic Timing Model 
(Two Post-Keplerian Parameters); Stellar Masses 

We undertook a series of timing analyses in which 
general relativity was assumed to hold, so that there 
were only two independent post-Keplerian parameters 
instead of the three parameters of the basic timing model. 
This timing model had a total of 22 parameters and 
5013 degrees of freedom. We chose the independent 
variables to be cosi and m2. We analyzed a uniform 
grid of values of these quantities, restricting cos i to be 
between and 1, and m2 to be between and 0.5 M©. 
For each value of cosi and m2, the values of Pb, mi, 
r, and s were calculated using equations through |S1 
The TOAs were fit to a timing solution with these 
parameters held fixed but all other parameters allowed to 
vary. For each combination of cos i and m2 , the statistic 



Ax' 



X 



Xn 



of the best fit timing solution was 



recorded, where is the goodness-of-fit of a particular 
timing solution and Xmin the global minimum from all 
grid points. Fits with acceptable values of Ax^ at the 
68% and 95% confidence levels are shown in Figures (Sb 
andU 

The values of cos z, mi, and m2 and their uncertainties, 
were calculat ed followed the pr ocedure outlined in 
Appendix A of lSolaver et al.l l|2002l ^. Essentially this was 
a Bayesian analysis with uniform priors in cos i and m2 . 
A probability was assigned to each grid point based on its 
Ax^. After suitable normalization, the probabilities of 
all points associated with a given range of cosi (or mi or 
m2) were summed to calculate a probability distribution 
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Fig. 3. — Constraints on cosi and m2. Dashed lines indicate values of mi according to equation |H] (a) Left plot: Constraints from 
the basic timing model, with three post-Keplerian parameters (orbital decay and two Shapiro delay parameters), cast into inclination and 
mass values via equations |21 through |S] (b) Right plot: Constraints from the general relativistic timing model, with two post-Keplerian 
parameters (inclination and secondary mass). In each plot, inner and outer contours correspond to 68% and 95% confidence limits. 



function. The confidence intervals derived from these 
probability distributions were: 

f 0.4ll° J^ (68% confidence) 
cosi = < , „■ _ (7) 
[ 0.41+°;2^ (95% confidence), 

J 2.1±O.2M0 (68% confidence) 
" I 2.ll°:^ M© (95% confidence), 

and 

J O.191±O.O15M0 (68% confidence) 
" I 0.19lt°:°33 Mq (95% confidence). 

3.4. Lower Limit on Pulsar Mass from Pf, Alone 

While there is no reason to doubt the detection of 
Shapiro delay, it is worth noting that the orbit decay 
measurement alone provides some evidence for a very 
massive neutron star. We calculated a probability 
distribution function for mi assuming a uniform a 
priori distribution for cosz, appropriate for randomly 
oriented orbits, and assuming P{, is drawn from 
the Gaussian distribution implied by its measurement 
uncertainty. The distribution of mi calculated under 
these assumptions gives lower limits mi > 1.75 M0 (68% 
confidence) and mi > 0.88 (95% confidence). 

3.5. Solar Wind 

Radio signals are dispersed by electrons in the 
interstellar medium and in the solar wind. Because the 
solar wind is variable and unpredictable, its contribution 
can degrade the accuracy of pulsar timing models. 
The problem is particularly acute for PSR J0751-I-1807 
because of its low ecliptic latitude. Dispersion by the 
solar wind imposed delays of up to ^12 fis in the TO As 
in our data set. We excluded all observations for which 
the line of sight to the pulsar passed within 15° of the 



sun. We found that changing this cutoff to 30° or to 0° 
had little impact on our results. 

In principle, observations at multiple radio frequencies 
at every epoch would allow the dispersion to be measured 
and removed; in practice, it is difficult to measure high 
precision TO As at two or more widely spaced frequencies, 
so it becomes necessary to use an analytic model of the 
solar wind and average the results over many epochs. 
We modeled the electron density in the solar wind as 
ne(r) = rto(lAU/r')^, where r is the distance to the 
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Fig. 4. — Constraints on pulsar and secondary masses from 
the general relativistic timing model. Confidence limits of 68% 
and 95% are shown. These are the same constraints as the 
right plot of figure |21 cast into a different parameterization. The 
shaded region in the lower left is disallowed by the Keplerian mass 
function. Dashed lines show constraints from alone. A dotted 
line indicates inclination i = 60° . 
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sun and uq is the electron density at 1 AU. We found 
the best fits to the data had no = 9.6 ± S.Opcem"^. 
Within this uncertainty range, the particular value of 
no used had little impact in the rotation and binary 
parameters derived in the timing fit; we used a fixed value 
of no = 9.6 cm"'^ to calculate the rotation and binary 
parameters in Table |21 

The pulsar's astrometric parameters — position, proper 
motion, and parallax — exhibit high covariances with 
no and with time changes in no. To estimate values 
and uncertainties for these parameters, we fit a grid of 
timing solutions with no between 6.6 and 12.6pccm~^, 
and with the time derivative of no between —1.5 and 
+1.5pccm^'^ yr^^. We use the extreme values from these 
timing fits to calculate the uncertainties on position, 
proper motion, and parallax given in Table [3 

Observations of another pulsar, J1713+0747, over a 
similar period of time found a marginally smaller solar 
wind electron density, no = 5 ± 4 pc cm~'^ ijSplaver et al.l 
l200f^. Despite differing by a factor of two, the no values 
from the two pulsars are in statistical agreement. 

4. PROBING STRONG FIELD GRAVITY 

Pulsars are well established testbeds for relativity. 
Observations of the neutron star-neutron star binary 
PSR B1913+16 have established that its orbit decays 
at the rate predicte d by g ei ieral relativity with in 0.3% 
lITavlor fc Weisberd IT9M IWeisberg fc Tavlod l20(il . 
However, this test of relativity is of limited use for 
constraining violations of the strong equivalence principle 
(SEP) because of the similar self-energies of the two 
neutron stars. More useful for probing this aspect of 
gravitation are neutron star-white dwarf binaries, in 
which the radically different self-energies of the two 
stars would generate excess gravitational wave energy 
loss under SEP-violating theories. A succinct review 



is given by 'Arzoumanian (2003) fsc 
Will & Zaelauer.. l989; Goldman.,1992 
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The observed change in the orbital period 
PSRJ0751-hl807is iPb/Pb)ohs = (-2.7±0.4)xl0-i«s 
We consider five possible mechanisms for generating 
changes in the orbital period: acceleration of the binary 
system relative to the Earth; changes in the gravitational 
constant, i.e., nonzero G; changes in the masses of 
the component stars; energy loss from the binary due 
to quadrupole gravitational radiation, as in general 
relativity; and energy loss due to dipole radiation, as 
in SEP-violating theories. Contributions from tidal 
interactions between the two stars arc likely not to be 
important (e.g., Smarr & Blandford 1976). The observed 
decay rate is the sum of the five contributions: 



(10) 



We address each of these terms in turn. 

The first term, acceleration biases, arises from the 
proper motion of the pulsar, acceleration toward the 
Gal actic plane (z-accelerati on) , and Galactic rotation 
(see|D amour fc TavlodlTggil for details). The biases for 
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PSRJG751+1807 are listed in Tabled These values 
were calculated using the measured proper motion of 
6masyr~^, and using a distance of l.lSkpc, derived from 
the dispersion measure of the pulsar using the NE2GG1 
galactic electron density model (Cordcs & Lazio 2002). 
To be conservative, we assign an uncertainty equal to 
the total bias (i.e., 100% uncertainty). The bias of Pt 
is much smaller than its measurement uncertainty and 
essentially negligible.^ 

The second term in Equation ^1 due to G, can 
be shown negligible by appeal to other binary pulsars. 
The change in orbital period due to nonzero G is 
approximately {Pb/Pb)Q — —2{G/G). This expression 

neglects the effects of G on the energy content of the 
stars themselves; a more precise expression is given in 
iNordtvcdt (1990), but is not necessary for the purposes 
of the present paper. Because the same expression 
for (Pb/Pf,)(5 holds for any binary pulsar system, the 
binary with the lowest value of this expression sets 
an upper limit on it for all pulsars. The best limit 
comes from PSR J1713+0737, for which Pb/Pb < 1 x 
lO-i^s-i (ISnlaver et a] . 1 1200,511 . more than an order of 
magnitude smaller than our value for PSR J0751+1807. 
A comparable limit on (Pb/Pb) q can be set using limits 
on G from lunar laser ranging measurements, which have 
found G/G = (1.3 ± 2.9) x lO-^Os"! llWilliams et alJ 
120041) . In any case, G effects are unimportant for 



PSRJ0751+1807. 

The third term arises if one of the component stars is 
losing mass ijEsposito fc Harrisonlll975t iLavaeetto et alJ 
l20?)M . First, we consider mass loss from the pulsar. 
The pulsar's measured spin down rate implies that it 
is losing energy at a rate E — Att^vqViIi, where Ii is the 
moment of inertia of the pulsar. Assuming the pulsar 
outflow is relativistic, this energy loss rate is equivalent 
to a mass loss rate of nii = E/c^ = Att^voVxIx/ = 
-4 x 10-21 MQS-i(/i/10'*'^gcm). The resulting change 

in the orbital period is {Pb/Pb)m = — 2ni,i(mi +m2)^^ ~ 
—4 x 10-21 s-i, depending on the precise values of nii, 
7712, and Ii. This is three orders of magnitude smaller 
than the measured value and, so mass loss can be 

^ Tlie same phenomena shift the pulse period derivative, P = 
ui/uq, a few percent away from its intrinsic value, but this is of no 
practical consequence. 
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neglected. 

Next, we consider mass loss from the white dwarf 
companion. I. Wasserman (private communication) has 
pointed out that mass flow from the white dwarf could 
arise from irradiation of the white dwarf by the pulsar. 
If such an outflow were directed predominantly opposite 
the motion of the white dwarf, it would induce a negative 
Pfc. Given the geometry of the PSR J0751+1807 system, 
a white dwarf of low mass might capture sufficient 
the pulsar flux to produce such an outflow, although 
there is no reason to expect the flow to be collimated 
tangential to the orbit. The presence of an ionized 
outflow would give rise to variations in the dispersion 
measure of the pulsar over the course of its orbit and, 
depending on the orbital geometry, might even cause 
eclipses. However, PSR J075H-1807 exhibits neither 
eclipses nor orbital variability in its dispersion measure. 
To search for the latter, we analyzed multi-frequency 
observations made in January 2004 by dividing the orbit 
into ten equal length sections and separately measuring 
the dispersion measure in each section. We found 
no significant variation in dispersion measure over the 
course of the orbit, with a conservative upper bound of 
ADM < 4 X 10-4pccm-3 1 x lO^^cm-^. This is 
strong, but not definitive, evidence against an outflow; 
it is possible to imagine scenarios with small orbital 
inclinations in which the pulsar irradiation directs the 
outflow away from our line of sight. Optical observations 
of the companion find it to be relatively cool, further 
evidence that a significant wind is unlikely (van Kerkwijk 
et al. 2005; M. van Kerkwijk, private communication). 

The expressions for {Pb/Pb)Q and {Pb/Pb)D in 
Equation 1101 depend on gravitational theory. In general 
relativity, {Pb/Pb)Q is given by Equation |21 while 
{Pb/Pb)D ~ 0. In some SEP-violating theories, 
the quadrupole term equals the general relativistic 
expression, while the dipole contribution has the form 

\Pb J ^ \PbJ \GJmi+m2 

(11) 

where ad and ac2 are the couplings of the pulsar and 
the secondary star, respectively, to the scalar field, and 
G^, is the "bare" gravitational constant. 

To measure or constrain the dipole term {Pb/Pb)D 
using Eguation llOl it is necessary to know the quadrupole 
term {Pb/Pb)Q, which in turn requires knowing mi 
and 7712 by means other than the measured Pb. 
For PSR J0751+1807, there are no independent high 
precision measurements of 777i and 7712, so we must 
consider all combinations of mi and m2 which fall within 
the Shapiro delay 95% confidence contour (Fig|3Ji) and 
have reasonable neutron star m asses, 1 Mq < mi < 
SM© (|Lattimer fc Prakas'hll2004|) . For these masses, the 
predicted general relativistic values of {Pb/Pb)Q range 
from -1.1 X 10^18 s-i to -1.4 X lO^^'^s^^ The largest 
negative allowed value of {Pb/Pb)D — {Pb/Pb)ohs — 
{Pb/Pb)Q, is -4.1 X lO-^^s"^ (95% confidence). This 
arises at the smallest values of mi and m2 that fall 
within the constraints, and corresponds to a difference 
in coupling strengths 

(aci-ac2)'<7xl0-^ (12) 



This improves by a factor of a few on previously 
pubhshed up per limits fro m PSRs B0655-I-64 and 
J1012+5307 lArzoumanial I200I iLange et alJ l200l|) . 
A result comparable to it can be derived from the 
measuremen ts of the PSR J1141— 6545 orbit presented in 
iBailes et al.l l|2003 1 . In fact, in generalized tensor-scalar 
theories of gravity, the latter pulsar is likely to be 
the most constraining among known pulsar-white dwarf 
binaries. Its eccentric orbit gives rise to several 
relativistic phenomena, allowing its stellar masses to be 
measured independently of Pb ((Esposito-Farese .200^ . 

In Brans-Dicke theory, the coupling to star i is 
aci = 2si(2 -t- wbd)"^, where Si = —dlnmi/dlnG, the 
"sensitivity" of star i, is the change in its binding energy 
as a function of G, and where wbd is the Brans-Dicke 
coupling constant. Estimates of si range from 0.1 to 
0.3, depending on the neutron star equation of state 
(|Will fc Zaglaueilll989() . The sensitivity of the white 
dwarf is negligible. The limit on coupling strengths 
(Eqn. [m places a lower limit wbd > 1300(si/0.2)2 
on the Brans-Dicke coupling constant. This limit is 
higher than previous constraints from pulsar work, but 
lower than the best limit attained by other means, 
ynn > 40000, from radio ranging to the Cassini probe 
lIBertotti et a1jl200l. 

5. DISCUSSION 

The mass of PSR J0751+1807 is the largest measured 
for any pulsar. As shown in Figure [SJ pulsars in circular 
orbits with helium white dwarf companions tend to have 
masses greater than the canonical value of 1.35 M©. 
This is in contrast to pulsars and secondary stars in 
neutron star-neutron binaries, which fall within the 
range 1.18— 1.44M0 The relatively high masses of 

pulsars in neutron star-white dwarf systems presumably 
result from extended mass accretion during the late 
stages of their evolution. 

An inverse correlation between orbital period and 
pulsar mass is apparent in Figure [3 Any such 
relation is likely to be complicated by the different 
evolutionary paths followed by different systems. The 
four systems with the longest orbital periods, those 
of PSRs J1713-h0747, B1855-H09, J0437-4715, and 
J 1909— 3744, are classical wide millisecond pulsar-helium 
white dwarf binaries, which underwent extended stable 
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Fig. 5. — Measured pulsar masses in circular pulsar-helium 
white dwarf binary systems as a function of orbital period. Data 
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Ivan Straten et all 1200 li) . and .Nice et al. C2005) . and references 
therein. 
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mass transfer . In the library o f evolutionary tracks 
calculated bv iPodsiadlowski et alJ (^002), such systems 
show an inverse correlation between orbital period and 
pulsar mass, although there is a wide spread in pulsar 
mass for a given orbital period. In these calculations, it 
was assumed that half of the mass lost by the secondary 
was accreted onto the neutron star; whether the mass 
accreted onto the neutron star is, in fact, proportional to 
the mass lost by the secondary remains an open question. 

PSRs J0751+1807 and J1012+5307, with their 
short orbital periods, provide a challenge to binary 
evolution theories. Evolution calculations tend not 
to produce systems with binary periods close to the 
orbital period of PSR J0751+1807 at the end of mass 
transfer. Its orbital period falls above the periods of 
ultracompact systems but below th e periods of standard 
pulsa r- white dwarf b i naries (e.g., IPodsiadlowski et aP 
|2002() . lErgma et all l)2001|) are able to produce the 
properties of the PSR J075H-1807 system by invoking 
magnetic braking and heating of the secondary by 
irradiation from the pulsar. The latter effect drives mass 
loss in the secondary, increasing the orbital separation 
and preventing the binary from shrinking into an 
ultracompact system. This picture is supported by 
optical observations, which find that the secondary lacks 
a hydrogen envelope and that it has cooled rapidly 
l|van Kerkwiik et alJl20?)l . 

6. SUMMARY 

The orbit of PSR J0751+1807 decays at a rate of A = 
-(6.4 ± 0.9) X 10"". This is in line with the expected 
value from general relativity. Combined with Shapiro 
delay measurements, it implies the pulsar and secondary 
star masses are 2.1 ± 0.2 Mq and 0.191 ± 0.15 Mq, 
respectively. The mass of PSR J0751+1807 is the largest 
recorded for a pulsar, and it may imply greater mass 



is transferred in tighter low mass neutron star binary 
systems than in wider systems. 

The maximum mass attainable by a neutron star 
depends on the stars comp osition and the equa t ion o f 
state of nuclear matter. iLattimer fc P rakash (2001|) 
calculated neutron star mass-radius relationships for a 
number of plausible equations of state. They found 
that models with exotic components, such as pure quark 
stars or mixed phases with kaon condensate or strange 
quark matter, allow neutron stars to attain masses no 
higher than ~ 2M0. While the measurement of the 
PSR J075H-1807 mass is nominally above this limit, the 
measurement uncertainty is not yet small enough to draw 
firm conclusions on this point. 

The uncertainty in the measurement of Pf, scales 
with observation span to the —2.5 power for uniformly 
sampled data. The precision of this measurement can, 
therefore, be increased dramatically with a relatively 
small number of observations over the next few years. 
The Shapiro delay measurement, because it has no such 
time dependence, will improve more slowly, and hence 
uncertainty in the Shapiro delay (and hence in the 
inclination) will ultimately dominate the uncertainty in 
the mass measurement. 
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